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The neutron pairing correlation and the soft dipole excitation in medium-mass nuclei near drip- 
line are investigated from a viewpoint of the di-neutron correlation. Numerical analyses by means of 
the coordinate-space HFB and the continuum QRPA methods are performed for even-even ^*~^'*0, 
50-58Qg^ and *°~*®Ni. A clear signature of the di-neutron correlation is found in the IfFB ground 
state; two neutrons are correlated at short relative distances ;$2 fm with large probability ~ 50%. 
The soft dipole excitation is influenced strongly by the neutron pairing correlation, and it accompa- 
nies a large transition density for pair motion of neutrons. This behavior originates from a coherent 
, superposition of two-quasiparticle configurations \lx {1 + 1)]l=i consisting of continuum states with 

■ high orbital angular momenta I reaching an order of / ~ 10. It raises a picture that the soft dipole 

' excitation under the influence of neutron pairing is characterized by motion of di-neutron in the 

, nuclear exterior against the remaining A — 2 subsystem. Sensitivity to the density dependence of 

effective pair force is discussed. 
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C\| ! I. INTRODUCTION 

^ , Excitations of neutron-rich nuclei near drip-line are subjects which are currently investigated extensively. It is 
' expected that exotic properties such as halo, skin or presence of weakly bound nucleons in neutron-rich nuclei cause 
\f-^ ', new features in excitations. An example is the soft dipole excitation in light halo nucleiP,|2,ll,0,IEIE0l3i typically 
' in ^^Li and ^^Be, where significant El strength is observed above the very low neutron threshold energy. This is in 
00 ^ contrast to the situation in stable nuclei, where most of the El strength concentrates in the high energy region of 
Op . giant resonances. Recently the soft dipole excitation is observed also in heavier systems up to neutron-rich oxygen 
^ ■ isotopes 18-22 These oxygen isotopes do not exhibit noticeable halo structure [llj, but on the other hand 

, they contain many valence neutrons. This suggests that the soft dipole excitation is not always inherent to one- or 
' two-neutron halo, but instead it may be a many-body phenomenon more generally seen in many neutron-rich nuclei 
I . near drip-line reaching medium and possibly heavier mass regions. 

^ ■ The degree of collectivity or the nature of correlations responsible for this excitation is one of the central issues 
^ ' which need to be clarified. Indeed different mechanisms have been proposed so far. One of the simplest mechanisms 
' producing the soft dipole excitation is the one which is associated with uncorrelated excitation of a weakly bound 
' neutron to continuum states. In this case, significant El strength emerges just above the threshold energy of neutron 
escaping without resorting to correlations nor collectivity, as it is sometimes called the threshold effect ITsj . This 
arises from a large spatial overlap between the extended wave function of a weakly bound single-particle orbit, which 
^ is occupied by a neutron in the ground state, and that of low energy continuum orbits, to which the neutron is excited. 
■ - 1 Major aspects of the soft dipole excitation observed in one- neutron halo systems, e.g. i^Be, fit rather well with this 
uncorrelated excitation picture^,!^. However, when more than one neutrons participate in the excitation, correlation 
between neutrons plays a role, and different mechanisms of the soft dipole excitation can be expected. In the case of 
the two-neutron halo nucleus ^^L i, the p airing correlation among halo neutrons plays a decisive role for the binding 
and formation of the halo [3, [iM llfl Ii3j lill • It is suggested that the two halo neutrons in the ground state display an 
attractive correlation in such a way that they are spatially localized with respect to their relative distance, in a range 
smaller than the size of the nuclear matter radius 0, 0|. The pairing correlation of this type or the di-neutron 
correlation in short is predicted to cause a strong enhancement on the soft dipole excitation in the two-neutron halo 
nucleus [T^ l20ll. or even to form a collective vibrational motion of the correlated halo neutrons against the rest 
of the system [l aj 21 j . Experimentally signatures of possible di-neutron correlation in the soft dipole excitation are 
obtained in iiLi[3|, but strong neutron-neutron correlations are not probed in other experiments[l|,|3]- 

The soft dipole excitation in medium mass nuclei has also been investigated theoretically, but with different view- 
points and results. The random phase approximation (RPA) calculations based on the Hartree-Fock models or the 
relativistic mean-field model predict that the soft dipole excitations has a character of uncorrelated neutron excita- 
tion carrying very little collectivity in the neutron-rich oxygen isotopes [2^ 123 . |23 | . Note however that these RPA 
models do not take into account the pair correlations among neutrons. The shell model calculation reproduces 
rather well with the experimental soft dipole strength in oxygen isotopes, but continuum effects are not included. 
Recently quasiparticle RPA (QRPA) calculations that include explicitly the neutron pairing corrections have been 
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Derformed |26ll27ll28l l29| . Most recent ones predict that the neutron pair correlation has a sizable effect on El strength 
of the soft dipole excitation in the medium mass region (28t.i29i |. This suggests that the neutron pair correlation is an 
important key to clarify the character of the soft dipole excitation in medium mass and heavier systems. 

In the present paper, we analyze in detail pair correlation effects in the medium mass region with Z = 8 — 28, taking 
proton (semi-)magic oxygen, calcium and nickel isotopes as representative examples. Motivated by the debate on the 
light two-neutron halo nuclei, we pay special attention to possibilities of the di-neutron correlation in these medium 
mass nuclei. We will conclude that some features of the di-neutron correlation are indeed present rather generally in 
the ground state of the medium mass nuclei. Furthermore the di-neutron correlation brings about a characteristic 
and strong influence also on the soft dipole excitation. 

Our analysis is based on the Hartree-Fock-Bogoliubov (HFB) method in the coordinate-space representation for 
description of the ground state, and the continuum quasi part icle random phase approximation (the continuum QRPA) 
for the excitations. The coordinate-space HFB theorv|30l IsiL Is^ is beneficial for description of nuclei with pair 
correlations near drip-line in medium and heavier mass regions since it extends the general HFB theoryjs^ in such a 
way that contributions of weakly-bound and continuum orbits to the pair correlation is precisely treated on the same 
footing with all other nuclcons through an explicit account of the coordinate-dependence of the pair potential and the 
quasiparticle wave functions. It is possible then to describe the excitations by generalizing the coordinate-space HFB 
to a time-dependent problem, i.e. by considering_a linear response of the system against the external perturbation. We 
have formulated recently a new QRPA method 34] along this line to describe excitations embedded in the continuum 
region above the threshold energies of nucleon escaping, which is an essential feature of nuclei near drip-line. This 
continuum QRPA has several distinctive features, compared to the continuum RPA0, js^j, which neglects the 
pair correlation, and to the conventional QRPA neglecting continuum effects. One of the most distinguished aspects 
is that it enables us to include correlations among continuum two-quasiparticle configurations corresponding to two- 
neutron escaping as well as to one-neutron escaping. As the model formulation is fully microscopic in treating all 
the nucleon degrees of freedom democratically, the di-neutron correlation, if present, emerges only as a consequence 
of the microscopic description. In order to probe the di-neutron behaviors in the ground state and in the soft dipole 
excitation we look into the two-body correlation density and the pair transition densities, which provide information 
on pair motion of neutrons. 

Numerical calculations are performed for even-even neutron-rich oxygen, calcium and nickel isotopes ^*^^'^0, 
50-58(^g^ and ^^"^^Ni near drip-line and for some more stable isotopes for comparison. In Section ^ we analyze 
the di-neutron correlation in the HFB ground state. In Section ITTll analysis of the soft dipole excitations using the 
continuum QRPA method is presented. Conclusions are drawn in Section Hvl We do not discuss in the present paper 
the low-lying dipole strength in stable nuclei with neutron excess, called pygmy dipole resonance [s^, IS |3^ ^ 
since the situations are different from those in near-drip-line nuclei on which we put focus in the present paper. 
Preliminary report of the present work is seen in Ref . |42j| . 

II. DI-NEUTRON CORRELATION IN THE GROUND STATE 

A. Coordinate-space HFB with density dependent interaction 

Wave functions of weakly bound neutrons in nuclei close to drip-line extend largely to the outside of nuclear surface 
due to the quantum mechanical penetration. Since associated neutron density is very low, the pair correlations in 
internal, surface, and external regions may be different. The coordinate-space Hartree-Fock-Bogoliubov theory'33,'3^ 
Is^ enables us to deal with this non-uniformity by utilizing an explicit coordinate representation for the Bogoliubov 
quasiparticles, which are the fundamental modes of the single- particle motion under the influence of pairing correlation. 
In this scheme the quasiparticles are able to have both particle and hole characters simultaneously, and accordingly 
they are expressed by the two-component wave functions 

where tr — ±i ^tjij&nd q = n,p represent spin and isospin. The quasiparticle states are determined by the HFB 
equation 

Hoq(j)iq{r(j) = E,q(j),q{r(T), (2) 

with 
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where Eig is the quasiparticle energy. The HFB selfconsistent mean-field Hamiltonian Hoq consists of not only the 
particle- hole part hg — Xq including the kinetic energy term, the Hartree-Fock mean-field and the Fermi energy Ag, 
but also the particle-particle part hg originating from the pair correlation. The mean-field Hamiltonian, hg and hg, 
are expressed in terms of the normal density matrix pq{ra,r' a') ~ ($o| <^')i'q{i'<^) |$o), the pair density matrix 
Pq{ra,r'a') = ($o| ipqi''''a')ipg{ra) |$o) = (-2cr) {^o\-ipq{r' - a')ipg{ra) |$o) and the effective nuclear force. We do 
not need the explicit form of the correlated HFB ground state |$o) since ground state expectation values of various 
physical quantities can be evaluated with use of the Wick's theorem for the quasiparticle annihilation and creation 
operators (3ig and f3j^ satisfying the vacuum condition Pig |<i>o) = 0, and with use of their relation to the nucleon 
annihilation and creation operators given by 

i 
i 

In the present work, we derive the particle-particle mean-field hg by using the selfconsistent HFB scheme. As the 
effective nuclear force responsible for the particle-particle part hg, called shortly the effective pairing force below, we 
adopt the density dependent delta interaction 0, Q 

Vpa^rir, r') = ]^V^{1 - P,) [l - 5{r - r'). (5) 

With this choice, the particle-particle part hg becomes a local pair potential Ag(r) — ^ ^ Pi^'''^ expressed 

with the diagonal pair density Pq{r) — Pq{i'(J,ra). The parameter po together with the total density p(r) — 
Pn{r) + pp{r) in Eq.Q controls the density dependence of the effective pairing force. There is no established knowledge 
on the density dependence of the pairing force, being under current investigations [isl li^ li^ . In the following analysis, 
we consider three cases: 1) the pairing with strong density dependence (or the surface pairing force in short), for 
which the parameter po is set to the central total density po = 0.19 fm~'^, 2) the density independent pairing force 
with the choice of 1/po = (the volume pairing force), and 3) the case of an intermediate density dependence with 
Po — 0.32 fm^'^ (the mixed pairing force). Since a recent analysis suggests that the mixed pairing reproduces better 
the odd-even mass difference in many isotopic chains [4^. we adopt the mixed pairing as a reference choice. The 
surface and volume pairing forces are employed to examine sensitivity to the density dependence. 

Concerning the particle-hole part hg of the HFB Hamiltonian, we replace it by a spherical Woods-Saxon potential 
model in the present work for the simplicity of numerical calculation. The use of the Woods-Saxon model makes 
the continuum QRPA calculation feasible. To our knowledge, there is so far no fully selfconsistent continuum QRPA 
calculation that utilizes the full HFB Hamiltonian derived from an effective interaction. The parameters of the Woods- 
Saxon model follow Ref. [s^ , which gives a reasonable description of the giant dipole excitation in doubly shell-closed 
stable nuclei, such as ^^O and ^°Ca. Neutron single-particle orbits near the Fermi energy are shown in Fig^ for 
representative isotopes. 

Assuming the spherical symmetry, we solve the HFB equation ^ in the radial coordinate for each partial wave. Here 
we adopt the radial mesh size Ar = 0.2 fm, and the box size rmax — 20 fm, following Ref. [33| . A cut-off quasiparticle 
energy E^aax = 50 MeV and the angular momentum cut-off Imax = 12 are used in summing up contributions of 
quasiparticle states when we evaluate the density matrices and associated quantities. Note that we include continuum 
quasiparticle states with E > |Aq| which lie above the threshold energy |Aq|. These continuum states, whose upper 
component ifi^igira) of the quasiparticle wave function extends to the outside of the nucleus without decaying, 
contribute to the correlated HFB ground state. To achieve the selfconsistency between the quasiparticle wave functions 
and the HFB mean-field Hamiltonian an iteration method (3(i| is used. 

The force strength parameter Vq in Eq.(|31) is fixed for each isotope chain so that the calculated neutron average gap 
(A„) = f Pn(r)An(r)dr/ J pn{r)dr |30l l34 l5ll| gives an overall agreement with the odd-even mass difference of the 
three-point formulajig, as shown in Fig|21 A common value Vq = —280 MeV fm'^ is adopted in the case of the mixed 
pairing. Note that we here improve the pairing force parameter, compared with our previous calculations [2^ Is^ 
where the conventional systematics A^yst = 12 /^/A MeV is fitted to determine Vq. The value of Agyst in the oxygen 
isotopes is larger by about 40% than the experimental odd-even mass difference (FiglSJ . In the following we shall use 
a larger value of Vq which corresponds to the value of A^yst in order to investigate dependence on the pairing force 
strength. 
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B. Two-body correlation density 



In order to analyze spatial behavior of the neutron pair correlation and to search for possible di-neutron aspects, 
we evaluate the two-body correlation density 

Pcorr,q{rcr,r'(T') = ($o| ^ r,)S{r' - rj)Sa,aSa^„' |$o> - Pq{ra)pq{r'a') (6) 

for the calculated ground state. This quantity displays correlation between two neutrons at positions r and r' with 
spins cr and a' . The spin anti-parallel (spin-singlet) configuration aa' =][ is responsible for the neutron pairing. 
In the above definition we subtract the uncorrelated contribution pq{ra)pq{r' a') in order to separate the change 
originating from the correlation. 

The two-body correlation density is also expressed as 

PcorrAr<y. r'o') = |p,(ra, r'~a')\^ - \pq(To, rV')P (7) 

in terms of off-diagonal parts of the pair and the normal density matrices. For the spin anti-parallel configuration, the 
first term of Eq.Q gives dominant contribution when the pair correlation is present. A different behavior is seen for 
the spin parallel configuration, where the two-body correlation density probes dominantly the Pauli repulsion effect, 
which is brought by the second term. In the following we concentrate on the spin anti-parallel neutron correlation. 
In displaying this quantity, we fix the position r' of one spin-down (a' =4) neutron (called the reference neutron 
hereafter) and plot it as a function of the position r of the other spin-up (cr =t) neutrons. Actual plots are made for 
the two-body correlation density Pcorr,n{T T7''"' V)l Pniy' i) divided by the neutron density p„(r' |) at the position of 
the reference neutron. This represents the conditional probability of finding neutrons at position r with spin a =| 
provided that the reference neutron is fixed at r' with spin a' =J,. This normalization removes the trivial radial 
dependence of the density that falls off exponentially as the reference neutron position moves to the exterior region. 
This facilitates comparison among different positions of the reference neutron. 

Examples of the two-body correlation density are displayed in Fig|31for nuclei near drip-line ^^O, ^*Ca and '^'^Ni. 
The reference neutron is placed at a position r' = (0, 0, z') along the z-axis, where z' is fixed the surface radius 
z' ~ Rsurf- We evaluate Rsurf by a position of the half central neutron density. To examine dependence on the 
reference neutron position, it is further displaced at an internal (z' = Rsurf — 2 fm) and an external (z' = Rsurf + 2 
fm) positions shifted by ±2 fm from the surface. The external position z' — Rsurf + 2 fm represents the neutron skin 
region, as the neutron density at this position is about ^ 1/30 of the central density (see FiglHJl in these isotopes. 
(Note that a criterion /3„(r)/p„(0) = 1/100 is sometimes adopted in the literature p9ll50l | to define the neutron skin 
thickness and to distinguish from a typical neutron halo, which emerges with a lower density p„(r)//9„(0) < 1/100.) 

It is seen that the correlation density Pcorr,n{i' i) exhibits a large and sharp peak in all the cases shown in 
FigEl The peak position almost coincides with the reference neutron position. The width £^d^ evaluated by the full 
width at half maximum, of the peak of the correlation density is only about 2 fm when the reference neutron is placed 
at the internal or the surface positions. This obviously means a strong concentration of the correlation density around 
the reference neutron as the width is smaller than the nuclear radius or the distance between the reference neutron 
and the nuclear center. In the region other than the first largest peak, the correlation density displays oscillatory 
behaviors, but its absolute value is much smaller than the first largest peak. The observed concentration of the 
two-body correlation density in the small region around the reference neutron indicates that two neutrons with spin 
anti-parallel (spin singlet) configuration has a large probability to come close at short relative distances |r — r'| < 2 
fm. It may be possible to regard this feature of the neutron pairing as that of the di-neutron correlation. 

It is possible to quantify the extent of the di-neutron correlation by evaluating the first largest peak in the two- 
body correlation density. Here we note that the first term in Eq. jT)) gives the dominant contribution to the correlation 
density. This term expressed in terms of the pair density matrix has a direct relation to the pair correlation. It 
is customary to regard the pair density matrix pn{ra,r' d')* — {'^Q\'ijj'l^{ra)ij}'l^{r' a') |$o) as the wave function of a 
neutron pair in the correlated HFB ground state. In this sense the first term |p„(rcr, r'a-')^ = p„(rcr, r'cr') with 
the opposite spins represents the probability distribution of a spin-singlet neutron pair. Using this quantity and 
normalizing, we can define a relative probability 

, , /|r-rl<,..P»(^t,r-'i)dr- 

for the spin-up neutron to exist within a distance from the spin-down reference neutron. The quantity p{rd) with 
a suitable value of measures the probability for the correlated neutron pair to form the di-neutron peak. We call 
it the di-neutron probability in the following. Calculated examples of the di-neutron probability are listed in Table 
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n We adopt = 2 fm for the internal and the surface cases, while for the external case = 3 fm is chosen to 
cover the large peak near the reference neutron (cf. Figs|31[7|and|HJl- The di-neutron probability amounts to 30-60%. 
As a reference, we compare with an estimate which would be obtained if the neutron pair probability p„(r |,r' |) 
distributed uniformly in the whole nuclear volume. This uniform limit, which we evaluate by replacing p„(r "f , r' |) 
with the neutron density p„(r), gives p{rd) = 0.20,0.10,0.06 for z' = R^urf - 2,+0,+2 fm in ^sq, 0.08,0.04,0.03 
in ^*Ca, and 0.05,0.03,0.02 in ^''Ni. In contrast to the uniform limit, the microscopically calculated values oi-pird) 
shown in Table exhibit an significant enhancement especially at the surface and the external positions. 

The di-neutron correlation emerges systematically at different positions inside and outside the nucleus, as seen 
in Fig|3| Inspecting in more detail, we find that, besides the width mentioned above, features of the di-neutron 
correlation vary with the reference neutron position. When the reference neutron is moved from the internal position 
(z' = i?s„r/~2 fm) to the surface [z' = Rsurf)^ the di-neutron correlation apparently enhances. This is also seen in the 
di-neutron probability (see the cases of ^^O, ^*Ca and ^''Ni in TableP), which shows an increase fromp(rd) « 30 — 40% 
at the internal position to ~ 50% at the surface position. Moving further toward the outside, as represented by the 
external position in the skin region z' = Rsurf + 2 fm, the concentration of correlation density around the largest peak 
is constantly quite large, keeping the di-neutron probability p{rd) ~ 50 — 60%. It is seen in this case that the peak 
position of the correlation density deviates slightly from the reference neutron position. The deviation is however 
within the di-neutron width ^d, keeping a large spatial overlap of the other neutron with the reference one. As the 
reference neutron moves far outside the nucleus, the deviation from the reference neutron increases further whereas 
the spatial correlation survives rather robustly even around z' ~ Rsurf + 3 to -1-4 fm in ^^Ca and ^''Ni. (An example 
of the correlation density with z' — Rsurf + 3 fm is shown in Fig01) This behavior also manifests itself in the pair 
probability p{rd)- In ^^Ca for example, this quantity with the reference neutron positions z' = Rsurf + 2,-1-3,-1-4 
and -|-5 fm reads p{rd) = 0.59,0.57,0.52 and 0.53 with = 3,4,5 and 6 fm chosen respectively to encompass the 
first largest peak whereas p{rd) with fixed rd(=3 fm) decreases as p{rd) = 0.59,0.38,0.18,0.07 at the same reference 
positions. We find a similar behavior in ^"^Ni while in ^^O the strong concentration of two-body correlation density is 
seen for z' Rsurf + 3 fm. The above observations indicate that the di-neutron correlation is most strong in the surface 
and the skin regions while the spatial correlation associated with the di-neutron behavior remains and decreases only 
gradually even outside the skin. 

If we compare the oxygen, the calcium and the nickel isotopes in Fig|21 we observe that the concentration of two- 
body correlation density around the position of reference neutron is more evident in ^^Ca and *'*Ni than in ^^O. For 
example, the di-neutron probability in ^^Ca and ^''Ni is as large as the one in ^^O (TableHJ. This means that the 
enhancement relative to the uniform limit is much larger in ^^Ca and especially in ^''Ni than in ^^O. It is also seen 
that the oscillatory behavior apart from the largest peak is weaker in ^^Ca and ^''Ni than in ^^O, especially when the 
reference neutron is placed at the external position. The small oscillation is a remnant of nodal structure in the wave 
functions of neutron single-quasiparticle states near the Fermi energy, e.g., 2si/2 and 1^5/2,3/2 neutron states in the 
case of ^^O although the magnitude of oscillation is suppressed by coherent contributions of other neutron quasiparticle 
states (see the next subsection for details). The smearing of the single-particle structures is more effective in heavier 
systems as more single-particle levels participate in the pairing correlation. It should be noted also that among the 
three isotopes the neutron separation energy (related to the Fermi energy) is smaller in ^^Ca and especially in ^^Ni 
than in ^^O (see Fig^. This difference in the neutron binding also influences the neutron pairing correlation in the 
external region as discussed just below. 

We have also analyzed the neutron two-body correlation density along the isotopic chains of Ca and Ni to check 
the di-neutron property in more stable nuclei and to examine how the di-neutron correlation varies with approaching 
the neutron drip-line. We select ^"^Ca and ^^Ni as examples representing stable nuclei and neutron-rich unstable 
nuclei (situated between the drip-line and the stable region), respectively. Here ^^Ni is chosen to compare with the 
near-drip-line nucleus ^*Ca having the same neutron number N — 38, and also with the near-drip-line isotope ^^Ni. It 
is found that the gross behavior of the two-body correlation density in ^''Ca and ®^Ni is similar to those in ^^Ca and 
^^Ni in the surface and the internal regions. As a representative example, we show in Fig|5lthe two-body correlation 
density in ''^Ca and ^^Ni for the reference neutron fixed at the surface {z' = Rsurf)- In Table^we do not see obvious 
difference between the near-drip-line nuclei and the more stable ones in the di-neutron probabilities p{rd) at the 
internal position z' — Rsurf — 2 fm and at the surface z' = Rsurf- 

A clear difference emerges, however, as the reference neutron is placed far outside the nuclear surface. This is 
illustrated by Fig01 where we compare the correlation densities in ^^Ni, ^^Ni and ^^Ca with a reference neutron 
fixed at a position (z' — Rsurf + 3 fm) in the far outside region. It is seen that the value of the correlation density 
evaluated at the reference neutron position is lower in ®^Ni by a factor of two or more than in near-drip-line nuclei 
^*Ca and *'*Ni. It is also seen that the large correlation density at the position of the reference neutron is most 
significant in ^'*Ni. This indicates that the di-neutron correlation in the external region is stronger in near-drip-line 
nuclei having shallower neutron Fermi energy than that in more stable nuclei with deeper Fermi energy (see Fig^for 
the single-particle energies and the Fermi energy in these nuclei). We can see directly this property also in the profiles 
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of the pair density Pni'r) of neutrons, shown in Fig|B| The pair density Pnif) has comparable or larger magnitude 
Pnir) ^ Pnjr) with the normal density p„(r-) in the skin region r ~ Rsurf + 2 fm and far outside. As discussed in 
Ref . |3ll| . the exponential tail constant of the neutron pair density /5„(r) in the asymptotic external region is related to 
the Fermi energy. The asymptotic tail develops as the drip line is approached, and is longer than that of the normal 
density Pnif)- Consequently the di- neutron correlation in the external low density region becomes relatively stronger 
in nuclei near neutron drip-line. 

C. Di-neutron correlation and single-particle configurations 

In obtaining the significant di-neutron correlation, it is essential to adopt a large model space of neutron quasiparticle 
states including continuum orbits up to large quasiparticle energy {Emax = 50 MeV in the present calculation) and 
large orbital angular momentum {Imax = 12). On the contrary, significant di-neutron correlation would not be 
obtained if we used a small single-particle model space. 

In the case of the oxygen isotopes, for example, one may consider that neutron single-particle orbits 2si/2 
1^5/2,3/2 lying around the neutron Fermi energy are those most responsible for the neutron pairing. (Note that the 
Woods-Saxon 1^3/2 orbit is not a bound orbit but a resonance close to the zero energy. See FigQ]) However, if one 
truncates to a single-j orbit or to orbits in one major shell, the di-neutron correlation never shows up. This is because 
the two-body correlation density would exhibit a mirror symmetry with respect to the x-y plane if only orbits with the 
same parity are taken into account. Apparently one needs more single-particle orbits with additional orbital angular 
momenta having different parities. 

In FigCI we show partial contributions of neutron quasiparticle orbits truncated with respect to the orbital angular 
momentum I of the orbits. Namely the two-body correlation density is evaluated by including the neutron quasiparticle 
states only up to a cut-off orbital angular momentum Icut- Inclusion of low angular momenta I = 0, 1, 2 covering at 
least all the neutron bound Woods-Saxon single-particle orbits in ^^O is insufficient to bring about the di-neutron 
correlation at the surface and the external positions; see the line corresponding to Icut — 2. If we add I — 3 fr /2,5/2 
orbits {Icut = 3), the di-neutron correlation becomes more visible, but an approximate convergence is achieved only 
by including further I — A — 5 orbits. It should be noted here that the neutron quasiparticle orbits with high angular 
momentum {I > 2) arc continuum states with £^ > |A„|. It is the pair correlation that makes these high-Z continuum 
orbits to contribute to the ground state pair correlation. Figdindicates that the contimmm high-/ orbits are important 
also for Ca and Ni isotopes, for which a larger value of angular momenta up to / ^ 6 — 8 at z' = Rsurf and 1^7—9 
at z' = Rsurf + 2 fm are necessary. 

The large contribution of the neutron high-Z quasiparticle orbits is related to the fact that the di-neutron correlation 
has a characteristic small width ~ 2 fm, which indicates an attraction within a short relative distance between two 
neutrons. It is illustrative to consider a pair of neutrons which are correlated in the relative s-wave at zero relative 
distance, and whose center of gravity is located off the nuclear center. A corresponding two particle wave function 
is written as ~ 5{r — r') oc X^hn ^m('")-^'™('"') where I is the orbital angular momentum about the nuclear center. 
This expression indicates that a coherent superposition of orbits with all orbital angular momenta I are necessary for 
such a correlation. On a similar basis, we can argue that a superposition of angular momenta up to Im is needed to 
describe the di-neutron correlation which is dominantly in the relative s-wave and has a width approximately given 
by id ^ '^T jlu {r being the radial position of the center of gravity of the pair). This estimate based on the observed 
di-neutron width ~ 2 fm gives a qualitative (though not precise) account of the maximum angular momentum seen 
in FigEI 

We have also examined contribution of the orbits by truncating the neutron quasiparticle states up to a cut-off 
quasiparticle energy Ecut- Results are shown in Fig|Sl The di-neutron correlation never appears only with the small 
model space including only up to Ecut = 5 MeV, which usually covers most of the quasiparticle states in one major 
shell. The results with cut-off of Ecut = 10, 20 MeV display only weak di-neutron correlation, and one needs neutron 
quasiparticle states at least up to ~ 30 MeV to obtain a qualitative account of the di-neutron correlation. As the 
neutron Fermi energy is small (— A„ — 3.54, 2.13, and 0.72 MeV in the case of ^^O, ^^Ca, and ^''Ni, respectively), the 
quasiparticle states that contribute to the di-neutron correlation are mostly those embedded in the continuum energy 
region E > |A„|. Most of the neutron continuum states are non-resonant states except for a few corresponding to 
the hole neutron orbits (e.g. lpi/2,3/2 and lsi/2) which have specific quasiparticle energies. The slow convergence 
with respect to the quasiparticle energy indicates that the non-resonant continuum states give non-negligible and 
accumulating contributions to the di-neutron correlation. It is noted that contribution of quasiparticle states with 
relatively small quasiparticle energy is important in the case of the external position. This may be related to the fact 
the pair correlation becomes weaker in the external region. 



7 



D. Dependence on pair interaction 

We have examined also dependence of the di-neutron correlation on the pairing force strength Vq. Here Vq is varied 
from the reference value Vq = —280 MeVfm^ (which produces (A„) = 1.5 MeV corresponding to the experimental gap 
in ^-^O) to an increased value Vq = —350 MeVfm'^ ii^n) = 2.8 MeV corresponding to the conventional systematics 
Asyst = 12/V^ MeV). Results are shown in Fig|51 It is noted that increase of the pairing force strength Vq enhances 
the di-neutron correlation, i.e., increases the intensity of only the first largest peak, rather than causing overall 
enhancement of the two-body correlation density. This indicates that the di-neutron behavior is indeed a correlation 
effect caused by the attractive interaction among neutrons. 

It is expected that the density dependence of the pairing interaction influences the di-neutron correlation since 
the pairing force of the surface and the mixed types give stronger neutron-neutron attraction in the surface and the 
external regions than in the interior. The volume pairing force (the density independent force) does not have this 
feature. To examine the influence of the density dependence, we perform calculations with use of the surface and 
the volume pairing forces. The pairing force strength Vq is chosen to reproduce the same average pairing gap (A„) 
calculated with the mixed pairing force. The calculated two-body correlation density is shown in Fig^| It indicates 
that the density dependence indeed affects the di-neutron correlation. Namely, in the case of surface pairing, the 
calculated two-body correlation density is significantly large in the surface (e.g. z' = Rsurf) and the exterior regions 
(z' = Rsurf + 2 fm) but it is weak in the internal region (z' = Rsurf — 2 fm). The volume pairing force does not exhibit 
such a position dependence. In the case of the mixed pairing force, the two-body correlation displays an intermediate 
feature as seen from comparison of Figs 171 and ITUl 

We emphasize also that the di-neutron correlation cannot be properly described if one adopts the schematic seniority 
pairing force used together with the conventional BCS approximation, where a constant pairing gap Ap is assumed 
instead of the selfconsistent pair potential A(r). This is illustrated in Fig llll where we present a BCS calculation 
obtained with use of the standard analytic expression of u, ^-factors, the Woods-Saxon single-particle energies, and 
the gap constant Aq. We include all bound and discretized continuum Woods-Saxon neutron orbits (obtained with 
the box radius Vmax = 20 fm), and use the same cut-off parameters Emax = 50 MeV and Imax = 12 as in the HFB 
calculation. The value of Aq is set to that of the average neutron gap (A„) obtained in the HFB calculation. It is 
seen that the BCS calculation causes a significant overestimate of the correlation density when the reference neutron 
is placed at external positions. In the case of surface and internal positions, the disagreement with the HFB is less 
serious, but profiles of the two-body correlation density are not well reproduced, as shown in Fig^J The unwanted 
overestimate in the external region arises from contributions of the discretized Woods-Saxon orbits in the positive- 
energy continuum region, for which the BCS approximation is known to cause unphysical correlation!^. If we neglect 
the discretized continuum Woods-Saxon orbits to avoid this difficulty, however, the BCS approximation produces a 
result (dotted line in Fig lll|l which is far off the HFB result, and the di-neutron correlation never shows up. 

E. Di-neutron size and coherence length 

The small width i^c; ^ 2 fm of the di-neutron correlation, or the di-neutron size in short, provides a characteristic 
length scale of the neutron pairing correlation. We note that this di-neutron size should be distinguished from the 
coherence length (the Pippard's coherence length), which plays a central role in the standard BCS theory of the 
metal superconductivity. The Pippard's coherence length for a superconducting uniform matter is given analytically 
as = Tivf/t^^ with vp and A being the Fermi velocity and the pairing gapjH^I- A direct application of to the 
nuclear case with typical values of vp and A gives the coherence length of an order of lOfm which is larger than the 
nuclear radius. If this estimate is adopted, the wave function of neutron pairs is meant to extend in the whole region 
of nuclear volume 53]. The present numerical investigation, however, reveals that the probability distribution of the 
neutron pair wave function is far from uniform, and a large probability ~ 50% is concentrated in a small region with 
short relative distances |r — r'| ^2 fm while the rest probability spreads in the whole nuclear volume. Our finding 
suggests that we may need at least two length scales (one is the di-neutron size and the other is the coherence length) 
to characterize the nuclear pairing correlation. 

We note also that the two-neutron correlation is calculated by Barranco et al. [l^ in a microscopic HFB description 
of the neutron pairing in the non-uniform low density neutron matter in a Wigner-Seiz cell with a immersed lattice 
'nucleus', approximating the situation of an inner crust of neutron stars. The authors point out that the root mean 
square relative distance weighted with the two-neutron probability density coincides approximately with the Pippard's 
coherence length ^p. The calculated two-neutron probability itself (Fig. 3 in Ref.|5J|), on the other hand, shows a 
distribution that forms a sharp and large peak at short relative distances |r — r'| ^ 2 fm, indicating a behavior similar 
to the di-neutron correlation discussed in the present investigation. This might suggest that the presence of the di- 
neutron correlation with a small length scale ~ 2 fm is a quite general feature of the nuclear pairing correlation. The 
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HFB calculation in Ref.'54'| adopts the finite range Gogny force as the efi'ective pairing force. Combining the results of 
Ref . (5 4] and ours, it can be suggested that the qualitative feature of the di- neutron correlation persists irrespective of 
detailed forms of the effective pairing force. The quantitative aspects however will depend on the effective interaction 
as we already discussed in the previous subsection. 



III. DI-NEUTRON CORRELATION IN THE SOFT DIPOLE EXCITATION 



A. Continuum QRPA description of the soft dipole excitation 

We first recapitulate briefiy the continuum QRPA methodjs^l which we adopt to describe the dipole excitation 
of nuclei near drip-line. It provides a fully microscopic description of a linear response of the nucleus excited by 
an external perturbing field with taking into account all nucleon degrees of freedom. It is formulated as the small 
amplitude limit of a time-dependent extension (TDHFB) of the coordinate-space HFB theory, which we utilize for the 
description of the ground state. Consequently the description is constructed in a selfconsistent manner. The linear 
responses in the normal and the pair densities are the basic quantities of the description. They are governed by the 
RPA density response equation called also the Bethe-Salpeter eauation[33l|: 

5~p+,,L{r,^) \ = dr' \ i?°^^(r,r',^) K:p,,,(r')<5p+.,L(r', c.)/r'2 . (9) 

Here the excitation with multipolarity L and frequency oj, and use of contact forces are assumed. The functions 
KphiT) and Hpair{f) represent the residual interaction associated with the density variations. An important feature 
of the present response equation is that we here include the particle-particle channel, i.e., the second and the third 
rows in Eq.(|5Jl containing 5p^^qL{r' ,uj) and (5/5_ ^/^(r', w), which correspond to the variations in the pair densities 
Sp±^q{rt) = 6 ^ J2a (V-'J('"o-)'0|(»"ct) ± ^pq{rd-)^pq{ra)'^ as well as the one in the normal density 5pq{rt) — 

5 {^{t)\'Y^^'4}\r(j)'4>{r) in the particle-hole channel, represented by the first row. To derive Eq.®, a linear 

perturbation in the time-evolving TDHFB state vector is considered. The products Hpair{r):pi5p±^qL{r,uj) 

of the residual interaction Upairif) and the pair density variations represent the dynamical change SA{rt) of the 
pair potential associated with the time-evolution. The two-point function i?Q^^(r, r', w) is the unperturbed response 
function for the three kinds of densities p{r) and p±{r), which are indexed by a and f3. Through the recursive relation 
for the density responses in Eq.®, the RPA correlations with infinite orders of the residual interactions are taken 
into account. The RPA correlation acting in the particle-particle channel, which is associated with the residual pair 
interaction Upair, may be called the dynamical pair correlaMon\2^. l29[ l84| . As the external field, whose radial form 
factor is represented by Vqf^ir), we consider the dipole operator 

7 AT 

= E(^^1m)(^0 - e- ^(ryi J(r,), (10) 

i(En i^p 

in which the spurious center of mass motion is explicitly removed. As the effective nuclear force, we adopt the density 
dependent delta forces. We employ the same pairing force Vpair used in the HFB description of the ground state 
to derive the particle-particle residual interaction Kpair- Thus the selfconsistency is achieved in treating correlations 
in the particle-particle channel. As the effective force responsible for the particle-hole correlation, we adopt a delta 
interaction of the Skyrme type 

Vphir, r') = (to(l + x^P,) + mi + x^P,)p{r)) 5{r - r'), (11) 

which gives a reasonable description of the giant resonances in closed-shell stable nuclei [s^- With this choice the 
residual interaction reads 

^lT:\r) = |(l_a:o) + |((5 + X3)p(r)-(2 + 4a;3)p,(r)), (12) 

>^t'\r) = io(l + |I) + ^(5 + X3)p(r), (13) 

^paAr) - y(i-^)- (14) 
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In the present continuum QRPA method, a special attention is paid for treatment of the continuum states that 
play essential roles for excitations embedded in the energy region above the threshold of nucleon escaping. To this 
end we evaluate the unperturbed response functions -Rq by means of an integral representation that uses a contour 
integral in the complex quasiparticle energy plane [33|: 

dEj2 ^^'^ll+i^^' {TrA^go^,,Ar,r',E + nu + ^e)Bpgo,,l,{r',r,E) 



-Tr^aSo.gy (r, r', E)Bpgo^gi'j, (/, r^E-huj- ie)} . (15) 



Here GQ.qij{E) = {E — Ho.qij)^^ is the HFB Green function in the partial wave Ij, which describes propagation of 
nucleons under influence of the pair potential A(r) and the particle- hole mean- field. We use the exact form of HFB 
Green function go,qij (E) js^ , which is given as a product of the regular and out-going solutions of the HFB equation 
(O so that QQ gij{E) satisfies the out-going boundary condition appropriate for continuum quasiparticle states. The 
combined use of the integral representation Ea. p5|) and the exact HFB Green function is the key ingredient of 
the present continuum QRPA method. Note that the present scheme takes precise account of two quasiparticle 
configurations where two nucleons occupy simultaneously continuum orbits since the two quasiparticles are both 
described by the out-going HFB Green function in Ea. (|15(l . Furthermore the particle-particle and the particle- hole 
correlations acting among such two-quasiparticle configurations are included through the density response equation jnj . 
The energy-weighted sum rule (the TRK sum rule) is satisfied within about one percent thanks to the selfconsistent 
treatment of the pair correlations and the use of the exact HFB Green function. Detailed derivation and other aspects 
of the continuum QRPA method are discussed in Ref . [33| . 

Some physical quantities are calculated directly from the solution of the density response equation. The El strength 
function for the dipole excitation is given by Sei{E = huj) — dB(El)/dE — — flmX^g/ '^'"'^gl;*i('')*'^Pgi=i(''j^)- We 
can also characterize the excitation mode by means of the transition densities. The particle- hole transition density 
for an excited state at the excitation energy Ei (i.e. with the frequency Ui = Ei/Ti) is given by 

Pfqir) = (<i>.|E^J(^^)V'.(^^)l*o)=r;M(r)<LW, (16) 



C 

PfgiW = -— ^Im(5pqL(r,w,;) (17) 



with use of the density response SpqL{r,uj) in the particle- hole channel. We normalize the transition density by a 
constant C so that the transition amplitude for the El operator, Miq = J drr'^v^^*{r)p^ll'j^{r), gives the El strength 

i3(El,0+s 1~) — J^^ SEi{E)dE — SMfy integrated over a certain energy interval around the excitation energy 
E ~ huji of the state under consideration. 

The present continuum QRPA also enables us to evaluate two kinds of transition amplitudes for pairs of nucleons: 

P^^'(r) = ($,| i,l{r ^)^lir j) |$o) = Y£,,ir)P!qW, (18) 
PtJ^ir) = ($,| ^,(r i)^,(r T) |<i>o) = Yl„j{r)P!ll{r). (19) 

These pair transition densities provide information of how nucleon pairs move in the excited state. The first pair 
transition density P[g{r) can be related to an amplitude to produce the excited state by adding (or transferring) at 
the position r a spin-singlet nucleon pair with the relative s-wave to the ground state of the A — 2 system. We call it 
the particle-particle or the particle-pair transition density hereafter. The second pair transition density P^^{r), which 
we call the hole-hole or the hole-pair transition density, on the other hand, is related to an amplitude of producing 
the excited state by removing a spin-singlet nucleon pair from the A + 2 ground state. These pair transition densities 
are calculated as 

Pl^Lir) = ^lm{Sp+,qL{r,u;,)-6p^^qL{r,u;,)), (20) 

Pl'qlir) = -^lm{6p+,qL{r,u;,)+5p^,qL{r,u;,)) (21) 

with use of the pair density responses Sp±^qL{r,uj) obtained in the density response equation. Here the same normal- 
ization constant C is adopted as for the particle-hole transition density. 

The force parameters of the particle-hole residual interaction is chosen as io — —1100 MeVfm'^, — 16000 MeVfm'^, 
xo = 0.5, and 2:3 = 1 taken from Ref. |35j| . To achieve an ap pro ximate selfconsistency in the particle- hole channel, 
a renormalization of the force strengths, i.e., to. 3 ~* / x to.slSSj is adopted so that the lowest energy dipole mode 
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corresponding to the spurious center of mass motion has the zero excitation energy. The same radial mesh as in solving 
the HFB equation Q is used to solve Eq.®. The maximum orbital angular momentum in the two-quasiparticle sum 
of Eq.Q is Imax = 12 in O and Ca isotopes, which is the same as in the HFB calculation for the ground state. In Ni 
isotopes we use a larger value Imax — 17 and a larger radius cut-off Vmax ~ 25 fm to achieve better convergence in 
the continuum QRPA calculations. In this case, the HFB calculation is performed with the enlarged Imax and rmax 
although we do not see any sizable influence in the ground state properties. The same mixed pairing force as that 
in the HFB calculation is adopted as the reference choice. We adopt a small imaginary constant with e = 0.2 MeV, 
which corresponds to smoothing of the strength function convoluted with a Lorentzian with FWHM=0.4MeV. The 
calculated results shown below differ from those in the previous analysis|j2^, where calculations are done for ^^O by 
using the surface-type pairing interaction (po = 0.16 fm~^) and a larger force strength Vq = —520 MeVfm'^ adjusted 
to the conventional systematics of the pairing gap A^yst = 12/vG4 MeV ('^ 2.5 — 2.8 MeV in the oxygen isotopes). 
Here the adopted value Vo = —280 MeVfm'^ with po — 0.32 fm~'^ is fixed so as to produce the odd-even mass difference 
(~ 0.5 — 2.0 MeV) as discussed in Section ^ Other numerical details are the same as in Ref . psf . 



B. El strength near neutron threshold energy 

The El strength functions Sei{E) = dB(El)/dE calculated for even-even isotopes 1^-240^ ^^•^'^•^^Ca, and so-s^Ni 
near drip-line are presented in Fig ll2l In all nuclides the strength function exhibits significant distribution of the 
El strength just above the threshold energy Eth,i = niin(_E,;„) + |A„| of one- neutron escaping and far below the 
giant dipole resonance energy [Egor ~ 20 MeV in O and ^ 15 MeV in Ni). This low-energy El strength increases 
significantly as the neutron drip-line is approached. In many isotopes the soft dipole excitation is situated also 
above the threshold energy Eth.2 = 2|An| of two-neutron escaping, which becomes low in nuclei near the drip-line 
and especially in the nickel isotopes. The strength distribution in the giant dipole resonance region shows also rather 
strong isotopic dependence in the case of the oxygen isotopes. The peak at the zero energy corresponds to the spurious 
center of mass motion, for which the energy weighted strength is negligible. In the following, we concentrate on the 
soft dipole excitations. 

The strength function in the region of the soft dipole excitation shows a smooth profile as a function of the excitation 
energy. It is not possible to evaluate the resonance width as the strength does not form a sharp resonance peak. The 
smooth profile implies that neutron escaping has a large influence on the soft dipole excitations. This is explicitly 
seen by comparing with a calculation neglecting neutron escaping, as is shown in FigElfor ^^O. This calculation is 
done with use of the discretized continuum quasiparticle states obtained with the box boundary condition (the box 
radius Vmax = 20 fm), instead of using the exact Green function with the out-going boundary condition adopted in 
the present continuum calculations. In the discretized calculation there are several discrete peaks in the energy region 
of the soft dipole excitation. The four major peaks seen in the interval = 5 — 9 MeV have different behaviors in 
the transition densities while the transition densities in the continuum calculation varies smoothly as a function of 
the excitation energy. The soft dipole excitation cannot be represented by one of these discrete peaks. 

The El strength in the soft dipole region is experimentally measured in neutron rich oxygen isotopes, and the energy 
weighted sum of the El strength below < 15 MeV is extractedl^. We evaluate the corresponding energy weighted 
sum from the calculated strength function, and list in Table IhI The experimental value of the energy weighted sum 
is about 8% of the Thomas-Reiche-Kuhn sum rule value in ^^'^^O, and about 12% in ^°0||9|. The calculation gives a 
fair agreement with the experimental data in ^*'^°0, but it overestimates in ^^O. The agreement may be improved 
by refining the model Hamiltonian, e.g., the Woods-Saxon parameters, or by adopting a selfconsistent Hartree-Fock 
potential which is expected to be more realistic. Note that the shell modelj25j and the phonon coupling model|2^ 
based on the Skyrme HF-I-BCS-I-QRPA approach exhibit a saturation around ^^O. 

Comparing among O, Ca, and Ni isotopic chains, we observe slightly different behaviors in the soft dipole excitation. 
In the case of Ni isotopes, the El strength is distributed at very low excitation energy as the neutron threshold energy 
is very low {Eth,i ~ 0.70 — 2.45 MeV). This can be related to the small neutron Fermi energy in the ^ = 80 — 86 nickel 
isotopes, which is only about —1.50 ~ —0.40 MeV. Note also that the single-particle energies of the most weakly 
bound neutron Woods-Saxon orbits 3s]^/2 and 2^5/2 in the vicinity of the Fermi energy are small; 635^^2,2^5/2 — 1 
MeV (Fignj. In Ca isotopes, the increase of the El strength above the neutron threshold energy Eth,i is not as steep 
as in O and Ni, but instead a small peak is formed at the energy which slightly deviates from the threshold energy. 
This is because there is no weakly bound neutron s orbit in the calcium isotopes while p orbits (2^1/2,3/2) participate 
instead (see the following section for details). 
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C. Pairing effects on dipole strength 

We analyze effects of the neutron pair correlation on the soft dipole excitation, which is the primary issue in the 
present investigations. To visualize the influence of neutron pairing correlation, we perform a calculation where the 
neutron pairing interaction is switched off, i.e. by setting Vq = 0. (For open sub-shell nuclei such as ^^'^oq very 
weak pairing interaction Vq = — 28 MeVfm"^ is used to guarantee a J'^ = 0+ configuration in the last j-shell orbit 
partially occupied in the ground state. This choice produces a so small average pairing gap (A„) < 0.1 MeV that the 
pair correlation effects are negligible.) Calculated results are shown by the dotted line in Fig^l We immediately see 
a sizable pairing effect on the El strength of the soft dipole excitation whereas there is essentially no effect on the 
giant dipole resonance. To focus on the soft dipole excitation, we show a magnified portion of the dipole strength 
function in Figll4l 

It is seen that the effect of neutron pairing correlation on the soft dipole strength varies depending on the isotopic 
chains and mass numbers, i) In the oxygen isotopes, the neutron pairing correlation increases significantly the El 
strength in the soft dipole region in 

^"'220 (and shghtly less clearly in ^^O). TableHIllists the S(E1) value of the soft 
dipole excitation integrated over an energy window of 4MeV width above the one neutron threshold energy Et^ i. 
The pairing correlation increases the -B(El) value by about 50-80 % in ^^~^^0. The energy weighted sum of i3(El) 
up to £^ = 15 MeV is also tabulated in Table ITU Influence of pairing on the energy weighted sum is relatively small 
in comparison with that on the i?(El) value. This is because the pairing effect becomes weaker at higher excitation 
energies. The pairing effect is negligible in the case of ^^O, but this is because the pair correlation itself is very small 
((A„) = 0.66 MeV) in this nucleus, ii) In the calcium isotopes, the pairing correlation does not enhance the magnitude 
of El strength. Instead it shifts the low- lying dipole strength up in the excitation energy and/or it suppresses the 
strength of the soft dipole excitation. The suppression in the El strength is about 20% in ^''Ca (see Table ^Oj. Hi) 
In the nickel isotopes, the neutron pair correlation either suppresses significantly the El strength at low energies (in 
^^Ni by about 30%, see Table|lll|) or modifies the shape of strength distribution (in ^°~^''Ni, see FigHH). After all, the 
pair correlation certainly influences the strength distribution of the soft dipole excitation, but it can either enhance 
(as in the case of 20,22q^ q^. g^ppj-ggg (^s in ^'*Ca,*^Ni) the El strength. 

To analyze the pair correlation effect, it is useful to decompose it into the static and the dynamical mechanisms |28l 
01 . Note that the pair correlation changes the static mean- fields in the HFB Hamiltonian, in particular the pair 
potential A(r). This mean-field effect modifies the ground state configuration and the single-particle excitation, 
through which the excitation properties are also affected. We call this mechanism the static pair correlation effect. 
On the other hand, the RPA correlation associated with the dynamical variation in the pair potential 5^{ruj) gives 
the additional pair correlation effect on the excitation, which we call the dynamical pair correlation effect. In other 
words, the dynamical pair correlation effect originates from the residual interaction taken into account in the RPA 
equation (|5J|, while the static effect are present even in the unperturbed response. To examine these pair correlation 
effects separately, we have performed calculations where the dynamical pair correlation is neglected while keeping the 
static pairing effects. Result is shown in Fig^jwith the dashed line. It is immediately seen that both the static and 
the dynamical correlations give considerable effects on the dipole strength in the soft excitation region. 

Let us first focus on the static pairing effect. It is seen that the static effect is a major part of the net pairing 
effect on the dipole strength, and produces qualitative trends of the strength function, although the dynamical effect 
cannot be neglected for quantitative description. To get more insight to the static effect, we look into the unperturbed 
strength function which is obtained by neglecting all the residual interactionsfFigll5|l. In the unperturbed strength, we 
separate contributions from different two-quasiparticle excitations by selecting a specified pair of angular momenta of 
the two-quasiparticle configurations. Taking ^^O as an example, where the pairing effect on the soft dipole excitation 
is large, we find that the static pairing effect increases the unperturbed El strength. This arises mainly from a 
contribution of a neutron two-quasiparticle configuration exciting the 2si/2 state and the continuum p states coupled 
to L = 1, abbreviated as [2s]^/2 '^P*]l=i or more shortly 2s]^/2 ^P* (the asterisk denotes the continuum states), as seen 
from comparison between the left top and the left bottom panels of Fig ll5l This quasiparticle excitation is available 
only by taking into account the ground state pair correlation. This is because the Woods-Saxon 2si/2 orbit located 
above the Fermi energy (Fig^ can be partially occupied only if the pair correlation is included. The contribution 
of the 2si/2 X p* configuration to the dipole strength is added with those of 1^5/2 x p* configurations, which are 
dominant ones when the pairing is neglected (see the lower panel of Fig ll5|) . In ^^"^"O, the increase of dipole strength 
due to the static pairing effect is similarly seen, but it not very large since the occupation of the 2si/2 orbit is smaller. 
Another aspect of the static pairing effect is that it pushes up the strength to slightly higher energy. This is because 
the pair correlation increases the energy of the two-quasiparticle excitation than that of the corresponding unpaired 
particle-hole excitation. 

The static pairing effect in Ca and Ni isotopes appears different from that in the oxygen isotopes. For example, 
it causes the large suppression of the El strength in ^^Ni, which is explained in terms of the neutron 3si/2 orbit. 
This orbit would be fully occupied by the last two neutrons if we neglected the pair correlation (cf. Fig^). Since the 
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3si/2 orbit is only weakly bound (the Woods-Saxon single-particle energy c^s^i^ = —0.78 MeV) and has an spatially 
extended wave function, the particle-hole excitations from this orbit to the continuum p orbits bring about a large 
El strength just above the one-neutron threshold energy, which is often referred to as the threshold strength (cf. the 
right lower panel of Fig ll5|l . Once the neutron pair correlation is included, the associated strength is reduced as the 
3si/2 orbit becomes partially occupied, making the contribution of the two-quasiparticle configurations 3si/2 x 
to the dipole strength significantly smaller. In Ca isotopes, the unperturbed strength in the soft dipole strength is 
dominated by contributions of the two-quasiparticle excitations 2p]^/2,3/2 x d* and 2px/2.3/2 ^ fFig ll5|l . The increase 
in the quasiparticle energy of the 2^1/2,3/2 state caused by the pair correlation pushes up the peak around E ^ 1 
MeV by about IMeV. 

We thus conclude that the static pairing effect shows variety in different nuclides, depending strongly on low-energy 
quasiparticle states (around the Fermi energy) which are quite sensitive to the pairing correlation. 

We then look into the dynamical •pairing ejfect. We immediately see that, in contrast to the static effect, the 
dynamical pairing effect has a systematic tendency to increase the dipole strength in all examples shown in Fig ll4l 
although the magnitude of the increase varies. The effect on B(E1) associated with the soft dipole excitation is shown 
in Tables ITIl and UTTl The largest effect on B(E1) amoimting to 10 - 40% is seen in ^^-^^O, and - 10 - 15% in 
the nickel isotopes. The increase of the strength due to the dynamical pairing effect is found in the previous QRPA 
calculations for oxygen isotopes [2^ 0, We here find that the increase due to the dynamical pairing correlation 
is universally seen in spherical nuclei near drip-line in the medium mass region. Note also that the enhancement due 
to the dynamical pairing effect has a similarity to the pair interaction effect predicted on the soft dipole excitation in 
the two-neutron halo nucleus ^^Li[l^. 

D. Transition densities: particle-particle dominance in the soft dipole excitation 

Characters of the soft dipole excitation can be clarified more directly by looking into the transition densities. 
Choosing a representative energy, we evaluate the particle-hole transition density pP^{r), the particle-pair and the 
hole-pair transition densities PPP{r) and P^^{r). They are plotted in Figs 1 1 61 [T7I and ITHl It is seen in the particle-hole 
transition density ff^{r) that in the external region {r ^ Rsurf) the soft dipole excitation has significant neutron 
amplitude whereas there is essentially no amplitude for protons, indicating that only neutrons are moving in the 
external region. At and slightly inside the surface, the particle-hole amplitudes of neutrons and protons have the 
same sign, but with the opposite phase to the external neutron amplitude, indicating that both neutrons and protons 
in this region move coherently against the external neutron motion. The external neutron motion is responsible for 
the soft dipole strength. This behavior of the particle-hole transition density in the soft dipole excitation is commonly 
seen also in other RPA and QRPA calculations without and with the pair correlations 1231 124. 1281 l29l| . 

A novel finding in the present analysis is that the neutron particle-pair transition density PPP{r) has very large 
amplitude in the external region, where the amplitude even exceeds that of the particle-hole transition density ff^ir). 
The hole-pair transition density P^^{r) on the other hand is the smallest among the three transition densities, and 
almost negligible in the external region. The relation |FPP(r)| > |p?'''(r)| > \P^\r)\ in the external region is seen m 
all isotopes in Figs 1 161 [T7I and IT51 except in and ^*Ca, where PPP{r) is still sizable. The dominance of PPP{r) 
indicates that the soft dipole excitation has a character of a particle-particle excitation: It is more appropriate to 
characterize the soft dipole excitation as a motion of spin-singlet neutron pair in the external region, than to describe 
as a simple particle-hole excitation of a neutron to continuum states. We also note that the particle-pair transition 
density PPP{r) displays a characteristic isotopic dependence. The large particle-pair amplitude in ^^'^''Ca and 18-22Q 
decreases with increasing the neutron number, while the particle-hole amplitude increases in the other way. A similar 
but slightly weak isotopic dependence is seen in the nickel isotopes. 

In order to reveal origin of the particle-particle dominance in the soft dipole excitation, we investigate influence of 
the pairing correlations on the transition densities. Comparing with calculations where the pairing effects are fully 
neglected (shown in Figs^| 1171 and [T51 bv the dotted line), we immediately see that the neutron pair correlation 
brings about the large particle-pair amplitude in the exterior region r > Rsurf- In FigstTCl . [T7I and IThI we show also 
results obtained by neglecting the dynamical pairing correlation. It is seen that both the static and the dynamical 
pair correlations are responsible for the particle-particle dominance. In particular, the dynamical pair correlation has 
a dramatic influence to enhance P^^(r) by a factor of about two or more in the external region. In the following we 
shall investigate in more detail the static and the dynamical effects separately. 

i) The static pairing effect: Generally low energy two-quasiparticle excitations, which are building blocks of low- lying 
excitations under the influence of pairing correlation, carry simultaneously particle-hole, particle-pair and hole-pair 
amplitudes. On the contrary, only the particle-hole amplitude would exist if the static pairing correlation were absent 
in the case of j-shell closed nuclei. As an example, let us consider the neutron two-quasiparticle excitations 2si/2 x P* 
and 1(^5/2 xp* (and 1(^5/2 x/*), which give dominant contributions to unperturbed strength function in ^^O (cf. Fig tT^ . 
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We first note that the two-quasiparticle configuration 2si/2 x p* would never participate in the soft dipole excitation 
if the pairing correlation were absent. With the pairing correlation included, the quasiparticle state 2s has a large 
amplitude both in the upper and lower components (corresponding to the particle and hole components, respectively) 
of the wave function since it is located near the Fermi energy. Accordingly the two-quasiparticle configuration 2si^2^P* 
brings a large amplitude both in the particle- hole transition density pP^{r) and in the particle-pair transition density 
PPP{r). (The hole-pair transition density P^^{r) is small since both 2si/2 and p* are located above the Fermi energy, 
and particle characters are dominant in these quasiparticle states.) The amplitudes PPP{r) and pP^{r) associated with 
the configuration 2si/2 xp* is especially large in the exterior region as the quasiparticle wave function of the 2si/2 state 
is spatially extended to the outside. The other dominant configuration 1^5/2 x p*{f*) contributes also to the particle- 
pair transition density in a similar way. The nodal pattern of the transition densities is consistent with the mixture of 
these two-quasiparticle configurations. The particle-particle character of these two-quasiparticle excitations in i^^^oq 
decreases as the Fermi energy (the neutron number) increases. In ^^O, the particle-particle character becomes small 
as the quasiparticle states 1^5/2 and 2si/2 both have a dominant hole character. The qualitative trends observed in 
the transition densities can be connected in this way to the properties of the relevant quasiparticle states. Similar 
mechanism are applied to the calcium isotopes, where the neutron 2^3/2,1/2 states play a central role. In the nickel 
isotopes, the relevant neutron quasiparticle states are 3si/2 and 2^5/2 (see FigstTland llSII . 

ii) The dynamical pairing effect: It is clear that the static pairing effects discussed above explains only qualitative 
aspects of the transition densities since the static pairing effect alone explains about a half of the particle-pair transition 
amplitude. The dynamical effect adds an essential enhancement in PPP{r) especially in the external region which is 
most relevant to the soft dipole excitation. The increase by a factor of 2 in the particle-pair transition amplitude 
ppp(^r) corresponds to an enhancement of a factor of about four in the strength of neutron pair transfer. Namely the 
characteristic particle-pair dominance of the soft excitation is strongly affected by the RPA correlations. This means 
that the calculated soft dipole excitation cannot be explained as a few representative two-quasiparticle configurations 
which we find responsible for the static pairing effect, such as 2si/2 ^P* or 1^5/2 'xp*{f*) in the case of O isotopes. This 
conclusion is quite different from that of the RPA calculations neglecting the pair correlations 22, 24], which predict 
the soft dipole excitation in neutron-rich oxygen isotopes as a non-collective independent particle-hole excitation of 
weakly bound neutrons. 

E. Di-neutron correlation in the soft dipole excitation 

We shall investigate nature of the large dynamical pairing effect on the particle-pair transition density PPP{r). 
The large amplitude in PPP{r) itself indicates that the pair correlation enhances the probability to find two neutrons 
participating in the soft dipole excitation at the same position r. It is then tempting to interpret it in connection with 
the di-neutron correlation which we found in the ground state. To check this viewpoint, we examine contributions of 
high-/ quasiparticle orbits to the soft dipole excitations. 

We have performed calculations where contribution of quasiparticlcs with high angular momenta is truncated in 
evaluating the RPA correlations. In practice, we put an upper cut-off Icut to the sum over orbital angular momenta 
I and V of the two quasiparticle configurations in the density response function, Ea. l|15|) . If the pairing correlation is 
completely neglected, the angular momentum of quasiparticles contributing to the dipole response is limited in a small 
range < / < /°„( = locc + 1, where locc is the largest orbital angular momenta of the occupied bound Woods-Saxon 
single-particle orbits. In the case of oxygen isotopes ^*~^^0, for instance, neutrons would occupy s,p and d bound 
orbits in the ground state in the null pairing case, and hence only the angular momentum combinations [s x p]l=i, 
[p X d]L=i and [d x /]l=i contribute. The cut-off /°„j = 3 is sufficient in this case. As the pairing correlation is taken 
into account, however, all combinations including [I x [l + 1)]l=i with / > are allowed to contribute. Note that 
these high-Z quasiparticle states are all continuum orbits. 

Results of the truncated calculations are shown in Fig ll9l It is seen that contributions from neutron high-/ quasi- 
particle states with / > locc are essential to produce the large enhancement in the particle-particle transition density. 
In ^^O and ^''Ca, the angular momenta up to / ~ 9 and / ~ 10 respectively are necessary to approach the final result. 
In ®^Ni, the orbits up to / 13 contribute in the external region up to r ^ 12 fm, but very high angular momenta 
/ > 13 still continue to influence in the far outside r ^ 15 fm. The particle-hole transition amplitude pP''-{r), on the 
other hand, is affected very little by the high-/ continuum configurations with / > locc, as seen in Fig ll9l 

We thus conclude that the neutron correlation responsible for the large enhancement of the particle-pair transition 
density PPP(r) in the soft dipole excitation is associated with a coherent superposition of a large number of neutron 
two-quasiparticle configurations with angular momentum coupling [/ x (/ -I- 1)]l=i involving up to large values of /. 
The accumulating high-/ contribution can be viewed as an evidence that two neutrons carrying the soft dipole mode 
are spatially correlated at short relative distance in such a way that we have seen the di-neutron correlation in the 
ground state (cf. Section^TOJ. This suggests that the soft dipole excitation is characterized rather strongly by motion 
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of a spin-singlet di-neutron in the nuclear exterior against the remaining A — 2 system. 

It is noted that much larger values of angular momentum contribute in the soft dipole excitation than in the ground 
state. In ^''Ni, for example, we need angular momentum up to ^ ~ 13 to achieve an approximate convergence around 
r = 7 — 12 fm where the pair transition density has the dominant distribution. This is because the di-correlation 
correlation in the soft dipole excitation takes place in much far outside of the nuclear surface than that in the case 
of the ground state. Note also that in '^''Ni convergence of high-/ contributions to the particle-pair transition density 
ppp(r) is slow in the very far exterior r > 15 fm even around the maximum angular momentum 13^/ < lmax{= 17). 
This is related to the fact that two neutrons can escape simultaneously in the nickel isotopes where the soft dipole 
mode lies above the two neutron threshold energy Eth.2- The oscillation of PPP{r) at r^lO fm indeed indicates a 
sizable two-neutron escaping. The slow convergence with respect to the angular momentum at very large distances 
suggests that there exists the pair correlation between the two escaping neutrons. 

F. Dependence on pair interaction 

Since the soft dipole excitation is strongly influenced by the neutron pairing correlation, sensitivity of the soft dipole 
excitation to the effective pairing force is expected. We shall examine this issue in connection with the di-neutron 
correlation. 

We have performed an analysis to see dependence on the force strength Vb for ^^O taken as an representative 
example. We increase the strength |Vb| of the mixed pairing force from the reference value Vq — —280 MeVfm'^ 
(corresponding to (A„) — 1.5 MeV fitted to the experimental odd-even mass difference) to Vq — —350 MeVfm'^ 
((A„) = 2.8 MeV representing the conventional systematics A^yst = MeV), as is done in Section ^ The 

calculated El strength function and the particle-pair transition density are compared in Figl^OI As expected, the 
increase of Vq modifies the strength distribution of the soft dipole excitation: The dipole strength calculated with 
Vo = —350 MeVfm'^ is larger than that with Vq = —280 McVfm"^. As seen in the right panel of the FiglSUI the increase 
is caused essentially by the dynamical pairing correlation. More dramatic sensitivity is recognized in the particle- 
pair transition density PPP(r), on which the dynamical pairing effect becomes significantly larger with increasing Vq. 
These observations are consistent with the picture that the di-neutron correlation manifests itself in the soft dipole 
excitation. It is also in accord with the finding in the previous section that the increase of Vq enhances the di-neutron 
correlation in the ground state (See. Section III DI and Fig|5J. 

We have also examined sensitivity to the density dependence of the pairing force. To this end, we performed 
calculations with use of the volume pairing (the density independent) force, and the surface pairing force with strong 
density dependence as is done in Sectional Calculated results are shown in Figsl21land l22l which may be compared 
also with those using the mixed pairing force, representing an intermediate density dependent force (cf. Figs ll4l [TBI 
1171 and [TH|l . Figsl21l and 1221 clearlv show that results are very different with different density dependences. It is seen 
that the dynamical enhancement of the particle-pair transition density PPP{r) in the external region (r > Rsurf) is 
much larger with the surface paring force than with the mixed and the volume pairing forces fFig l22|l . This trend 
is also seen in the El strength fFig l21l) . In the case of the volume pairing force, on the other hand, the dynamical 
pairing effect on the strength function becomes almost insignificant, and even the particle-pair transition density is 
not strongly enhanced. 

Note that the neutron- neutron attraction acting in the low-density external region is proportional to the value of Vq , 
which differs as |Vb| ~ 180,280,380 MeVfm^ for the volume, the mixed, and the surface pairing forces, respectively. 
Thus the above observation implies that the soft dipole excitation, especially the associated particle-pair transition 
density, is quite sensitive to the effective pairing force among neutrons moving in the low-density part outside the 
nuclear surface. This is of course related to the fact that the soft dipole excitation is essentially a mode carried by the 
correlated neutrons moving in the nuclear exterior. In addition, the strong sensitivity to the density dependence is in 
accord with a similar behavior of the di-neutron correlation in the ground state (cf. Fig llUI) . This again supports the 
picture that the soft dipole excitation is strongly influenced by the neutron pairing correlation of the di-neutron type. 

IV. CONCLUSIONS 

We have investigated the neutron pairing correlations and their influences on the soft dipole excitation in medium 
mass nuclei near neutron drip-line from the viewpoint of the di-neutron correlation. 

The analysis using the two-body correlation density has revealed the presence of spatial di-neutron correlation in 
the pair correlated ground state in nuclei near drip-line. It is found that correlated neutron pairs exhibit a strong 
concentration of the probability of about 30-60% at short relative distances | r — r' | ^ 2 fm, which is much smaller than 
the nuclear radius. This di-neutron correlation enhances in the surface and the skin regions in near-drip-line nuclei 
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although it also exists inside the nucleus and also in stable nuclei along the isotopic chain. The di-neutron correlation 
originates from coherent superposition of the single-particle (quasiparticle) orbits with large orbital angular momenta, 
which are embedded in the continuum energy region. 

We have analyzed the soft dipole excitation to search for the di-neutron correlation in this mode. It is found that 
the particle-pair transition density of neutrons in the soft dipole excitation is quite large outside the nuclear surface. 
This originates from the dynamical pairing correlation among neutrons moving in the external region, i.e., the RPA 
correlation for the excited state caused by the pairing interaction. Indeed the dynamical pair correlation is responsible 
for enhancing the particle-pair transition density by a factor of about two or more. This indicates that the soft dipole 
excitation under the influence of neutron pairing correlation has a dominant particle-particle character, rather than 
an uncorrelated particle-hole excitation from a weakly orbit to continuum orbits. We find also that two-quasiparticle 
configurations [I x {I + 1)]l=i involving continuum high-Z orbits up to around I ~ 10 accumulate coherently to bring 
about the large particle-pair transition density. This strongly suggests that there is the di-neutron correlation among 
neutrons participating in the soft dipole excitation. We are thus lead to a picture that in the soft dipole excitation 
a spin-singlet di-neutron move outside the nucleus against the A — 2 subsystem. Our analysis reveals also that 
the characteristic neutron pairing effects are sensitive to the density dependence of the effective pairing force. The 
influence of neutron pairing correlation on the dipole strength is sizable, but not always causes enhancement. The 
di-neutron correlation emerges more clearly in the particle-particle channel. 

We expect that the di-neutron correlation present in the soft dipole excitation may be most easily and directly 
probed in the two-particle correlation among two neutrons escaping from the excited state, or in transfer of neutrons 
to the excited state. These processes may also be used as a probe to study the density dependence of the nuclear 
pairing correlations. These issues are interesting subjects for future investigations. 
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TABLE I: The di-neutron probability p{rd) in ^^O, ^^^Ca and *''Ni near neutron drip-line, and in more stable ^^Ca. and ^^Ni. 
The reference neutron position r' = (0,0, z') is fixed at the surface { z' = Rsurf ), the internal {z' = Rsurf — 2 fm) and the 
external {z' — Rsurf + 2 fm) positions. The surface radius Rsurf defined by the half central density of neutrons is also listed. 
The di-neutron probability is evaluated with = 2 fm except in the external case where we use = 3 fm. See also the text. 
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11.8% 


19.8% 


B(E1) [e^fm^] {Eth.i <E< Eth.i 


+ 4 MeV) 




full pairing 0.188 


0.254 


0.393 


0.702 


no dynamical pairing 0.131 


0.205 


0.354 


0.694 


no pairing 0.104 


0.175 


0.235 


0.718 



TABLE II: The calculated energy weighted sum = dE' E' dB [El) /dE' of the El strength for the excitation energy below 

£■ = 15 MeV in the oxygen isotopes, and B(E1) = dE'dB (El) /dE'le^fm^] in an excitation energy interval of 4MeV above 

the one-neutron threshold energy Eth,i- The value of is given as a fraction to the Thomas-Reiche-Kuhn sum rule value 
^TRK- We list also the results obtained without the neutron dynamical pairing correlation, and those neglecting all the neutron 
pairing correlations. 



5«Ca 


54 Ca 


58Ca 


B2Ni 


B2Ni 


84n; 


««Ni 


B(E1) [e^fm^] 














full pairing 0.73 


1.48 


1.76 


1.39 


2.17 


3.01 


3.72 


no dynamical pairing 0.64 


1.40 


1.72 


1.15 


1.86 


2.69 


3.41 


no pairing 0.68 


1.86 


1.83 


1.27 


2.25 


3.21 


5.11 



TABLE III: The calculated El strength B(E1) = j ^^'^ dE' dB {El) / dE' [e^irc?] of the soft dipole excitation in Ca and Ni 
isotopes. In Ca isotopes the energy interval [Ei,E2\ is chosen with E\ = 5.9, 5.2, 4.8 MeV and E2 = E\ + A MeV to enclose the 
soft dipole peak \a A — 50, 54, 58, respectively. In Ni isotopes the interval with Ei = Eth,i and E2 = Eth,i + 5 MeV is used. 
See also the caption of Table HTl 
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FIG. 1: The Woods-Saxon single-particle levels of neutrons in ■^■^0,^'*Ca and ®^Ni. The single-particle energy is indicated in 
the parenthesis. The Fermi energy obtained by the HFB calculation is drawn with the thick dotted line. 
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FIG. 2: The calculated neutron average pairing gap (A„) in O, Ca and Ni isotopes, plotted with the crosses. The experimental 
odd-even mass difference evaluated with the three-point formula is also shown by the dashed line. The dotted line is the 
conventional systematics Asyst = 12/-\/34 MeV of the pairing gap. 
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FIG. 3: The neutron two-body correlation density pcorr,n{r t,r' i)/pn{r' [) in ^^O, ^**Ca and *'*Ni calculated with the mixed 
pairing force is drawn on the x-z plane. A contour plot of the same quantity is attached in the bottom of each panel where 
the interval of contour lines is 0.001 fm~^. The symbol "X" on the x-z plane indicates the position r' = (0, 0, z') of the 
reference neutron. The results for ^^O, ^**Ca and *^Ni are listed in the top, the middle and the bottom rows, respectively. In 
the middle column the reference neutron is fixed at the nuclear surface z' — Raurf whereas it is placed at the external position 
z' = Rsurf + 2.0 fm in the right column, and at the internal position z' = Rsurf — 2.0 fm in the left column. See Tabled] for 

the value of Rsurf- 




FIG. 4: The neutron two-body correlation density evaluated along the z-axis in near-drip-line nuclei ^*Ca and ^''Ni, and in a 
more stable isotope ^^Ni. The reference neutron is fixed at z' — Rsurf + 3 fm, which is indicated by the arrows. 



FIG. 5: The same as Figl^lbut for ^''Ca and for ^^Ni. The reference neutron is fixed at the surface position z' = Rsurf- 
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FIG. 6: The density Pnir) and the pair density /5„(r) of neutrons in ^'''^^Ca and ®®'®''Ni. The solid and the dashed lines 
represent Pn(r) whereas the dotted and the dot-dashed lines are for Pn{r). 
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FIG. 7: The dependence of the neutron two-body correlation density on the orbital angular momentum cut-off Icut = 0, 1 , 2, 3, • • • 
for ■^^O, "'^^Ca and **''Ni in the top, the middle and the bottom rows, respectively. The two-body correlation density is plotted 
along the z-ax.is. The full HFB result is drawn with the solid line. The reference neutron position is fixed at the internal 
{z' = Rsurf — 2 fm), the surface («' = Rsurf) and the external (z' = Rsurf + 2 fm) positions in the left, the middle and the 
right columns, respectively. The arrows indicate the reference neutron position. 
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FIG. 8: The dependence of the neutron two-body correlation density along the z-axis on the quasiparticle energy cut-off Ecut 
for ^^O, ''**Ca and *^''Ni, to which the top, the middle and the bottom rows correspond respectively. The reference neutron fixed 
at the internal {z' = Rsurf — 2 fm), the surface {z' = Rsurf) and the external {z' = Rsurf + 2 fm) positions in the left, the 
middle and the right columns, respectively. The solid line represents the full HFB result. 




FIG. 9: The dependence of the neutron two-body correlation density on the pairing force strength Vb of the mixed pairing force 
in ^^O. The solid line is the result with the reference value Vo — —280 MeV fm'^ corresponding to (A„) = 1.5 MeV whereas 
the dashed line is for Vo = —350 MeV fm^ corresponding to (A^) = 2.8 MeV. The reference neutron is placed at the surface 

Z — Rsurf • 
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FIG. 10: The neutron two-body eorrelation density along the z-axis in ^^O, ^**Ca and *''Ni calculated with use of the surface 
and the volume pairing forces having different density dependence. Here the dashed, the solid, and the dotted lines in each 
panel display this quantity for the internal, the surface, and the external positions of reference neutron. The pairing force 
parameters are: Vo = -375, -395, -385 MeVfm^ (in ^^O, ®®Ca and **Ni, respectively) and po = 0.19 fm"^ for the surface 
pairing, and Vo = -190, -178, -180 MeVfm^ (^^O, ^®Ca,**Ni) for the volume pairing. 
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FIG. 11: The neutron two-body correlation density along the a:-axis in the BCS approximation for ®*Ca with use of a constant 
pairing gap A„ = 1.27 MeV, plotted with the dashed line. The HFB result (the solid line) is also shown for comparison. 
Another BCS result using only the bound Woods-Saxon neutron orbits is also shown by the dotted line (labeled with BCS'). 
The reference neutron is placed at the external position z' = Rsurf -h 2 fm in the upper panel whereas it is at the surface 
z' = Rsurf in the lower panel. 
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FIG. 12: The El strength function in the neutron-rich even-even oxygen, calcium and nickel isotopes near drip-line, calculated 
with use of the mixed pairing force and plotted with the solid line. The result obtained without the pairing correlation (by use 
of a very weak pairing force Vb ~ 0) is also plotted with the dotted line. The arrows indicate the one- and the two- neutron 
threshold energies Eth,i and Eth,2- Note that Eth,i = Eth,2 in *^Ni and ^*Ni, where there are no bound quasiparticle states for 
neutrons. 
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FIG. 13: The El strength function in O, obtained by using the discretized continuum states and the box boundary condition, 
in the case of the mixed pairing force. See also the text. 




FIG. 14: The same as Fie ll2l but emphasizing the low excitation energy region and effects ol the static and the dynamical 
pairing correlations. For the latter purpose, the result obtained by neglecting the dynamical pairing correlation (while keeping 
the static pairing correlation) is plotted with the dashed line. 
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FIG. 15: The unperturbed El strength functions in O, Ca and Ni. The upper panels display the unperturbed strength 
functions obtained with the mixed pairing force whereas in the lower panels the results without the pairing correlation is shown. 
The partial strengths selected by a pair of the angular momentum quantum numbers of the two-quasiparticle excitations are 
also plotted. 
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FIG. 16: The particle-hole transition density pfg^('') (the top plot in each panel), the particle-pair transition density Pfqi^{r) 
(the middle plot) and the hole-pair transition density Piql {r) (the bottom plot) of neutrons for the soft dipole excitation in the 
oxygen isotopes, plotted with the sold lines. The particle-hole transition density p\q]^{r) of protons is also shown by the dashed 
line with wide intervals. For the pair transition densities P^^(r) and P^^(r), the neutron amplitudes calculated by neglecting 
the dynamical pairing effect are also displayed with the dashed line. The dotted lines represent the neutron transition densities 
calculated by neglecting all the pairing correlations. The arrow indicates the surface radius (the half density neutron radius) 
Rsurf - The selected excitation energy is £ = 10.0, 8.5, 7.0, 5.5 MeV for i*. 20,22, 24q ^ ^ r,^^ g y^^y ^j^^ ^-^^ 

no pairing calculation), which are indicated also in the figure. The B(£'l) value listed in Table UTI is used for the normalization. 
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FIG. 17: The same as Fig |16l but for the calcium isotopes. The transition densities are evaluated at the peak energy of the 
soft dipole excitation indicated in each panel. The El strength listed in Table ITTTl is used for the normalization. 
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FIG. 18: The same as Figs ll6l and 1171 but for the the nickel isotopes. The same excitation energy is used to calculate the 
transition densities with and without the pairing correlation. 
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r[fm] 

FIG. 19: The dependence of the particle-pair transition density r^P''^{r) (in the top panels) and the particle-hole transition 
density r^p''^'{r) (in the bottom panels) of neutrons on the cut-off orbital angular momentum /cut of the neutron quasiparticle 
states in '^'^O, -'^Ca and ^-^Ni. The results with = lcuti= 3), 4, 5, •••,9 for ^'^O, h^t = lcuti= 4), 6, 8, 10 for '^^Ca, and 
Icut = lcut{= 5), 7, 9, 11, 13 for ®'*Ni are shown by thin lines. Here the volume element is multiplied to magnify the amplitude 
in the external region. For reference sake, the result with the full pairing eflfects and the one without the dynamical pairing 
correlation are shown by the thick solid and the thick dashed lines, respectively. 
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FIG. 20: (a)Dependence of the El strength function on the pairing force strength Vo of the mixed pairing force in 2^0. The 
sohd and the dotted Hnes are the results for Vo = —350 MeV fm^ and —280 MeV fm^, corresponding to the average neutron 
pairing gap (A„) = 2.8 MeV and l.SMeV, respectively. In the right panel, the result obtained by neglecting the dynamical 
pairing effect is plotted with the dashed line to visualize the dynamical pairing effect in the case of Vb = —350 MeVfm^. (b) 
The same as (a), but for the particle-pair transition density P^^{r) of neutrons. 
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FIG. 21: The El strength functions in ^^O, '^^Ca and *''Ni calculated with use of the surface and the volume pairing forces 
having different density dependence (the solid line). The left panels show results with the surface pairing with po = 0.19 fm~'^ 
while the right panels are those with the volume pairing. The results with the full pairing effects is plotted with the solid line 
whereas the dashed line represents those without the dynamical pairing effect. See the caption of Fig llOl for the adopted force 
parameters. 
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FIG. 22: The same as Fig l21l but for the particle-pair transition density P^^{r) of neutrons associated with the soft dipole 
excitation at the excitation energy indicated in the figure. 
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